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ABSTRACT 

One of the major differences between various explosion scenarios of Type la supernovae (SNe la) is 
the remaining amount of unburned (C+O) material and its velocity distribution within the expanding 
ejecta. While oxygen absorption features are not uncommon in the spectra of SNe la before maximum 
light, the presence of strong carbon absorption has been reported only in a minority of objects, typically 
during the pre-maximum phase. The reported low frequency of carbon detections may be due to low 
signal-to-noise data, low abundance of unburned material, line blending between C II A6580 and Si II 
A6355, ejecta temperature differences, asymmetrical distribution effects, or a combination of these. 
However, a survey of published pre-maximum spectra reveals that more SNe la than previously thought 
may exhibit C II A6580 absorption features and relics of line blending near -6300 A. Here we present 
new SN la observations where spectroscopic signatures of C II A6580 are detected, and investigate the 
presence of C II A6580 in the optical spectra of 19 SNe la using the parameterized spectrum synthesis 
code, SYN0W. Most of the objects in our sample that exhibit C II A6580 absorption features are of 
the low-velocity gradient subtype. Our study indicates that the morphology of carbon-rich regions is 
consistent with either a spherical distribution or a hemispheric asymmetry, supporting the recent idea 
that SN la diversity may be a result of off-center ignition coupled with observer line-of-sight effects. 
Subject headings: supernovae: general supernovae: individual (SN 2010Y, 2010ai, PTFlOicb) 



1. INTRODUCTION 

The typical pre-maximum SN la spectrum consists of 
overlapping P-Cygni profiles of intermediate-mass ele- 
ments (IMEs) and Fe-peak elements (IPEs) that indi- 
cate expansion v elocities on the order of 10 4 km s~ 4 
(jFilippenkol I1997D . Photometric properties of the rise, 
peak, and decline of a SN la light curve can be ex- 
plained by assuming that a substantial amount of 56 Ni, 
synthesized in the explosion, powers th e SN luminosity 
(Colgate fc McKed 119691: lArnettl I1982D . Consequently, 
results of spectroscopic and photometric studies have 
supported the idea that SNe la are the outcome of 
a thermonuclear explosion of a C+O white dwarf in 
a binary system (IHovle fc Fowler] 119601: INomoto et all 
1984 lElias et all [19851: llbenl |198 8t INomoto et all 120031: 
Chen fc Lill2009l : lHowellll2010f) . 

The explosion mechanisms that have been proposed 
differ by how the thermonuclear flame is propagated 
through the star's interior; i.e., a sub-sonic deflagration 
via thermal conductivity or a super-sonic detonation due 
to strong shock burning. Pure detonation models ap- 
pear unlikely since they conflict with observations by 
producing too much 56 Ni, not enough IM Es, and leav- 
ing very little unburned material behind (jArnettl [l969t 
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IBranch fc" Khokhlovl I1995T) . On the other hand, while 
pure deflagration models may account for fainter SN la 
events, they are energetically weak, leaving too much 
material unburned to represent the majority of SNe la 
(iTrayaglio et al.H2004t IGamezo et alll2005l ; IKozma et all 
120051) . Thus, a deflagration that transitions into a deto- 
nation may be necessary in order to attain the nucleosyn- 
thetic yields that are consistent with the observations 
(IKhokhlovl[T99H IHoflich et all [19951: iKasen et al.ll200l 
iMaeda et alll2010al) . 

Additional facets of modeling (e.g., multidimensional 
consi derations) will affect the abundance tomography as 
well (|Wooslev et al.ll2009D . This can make it difficult 
to distinguish the dominant sources of SN la diversity. 
For certain, however, one similarity between all explosion 
models is the existence of burned (ash) and unburned 
material (fuel). By using spectroscopic signatures of C 
and O at pre-maximum epochs, one can infer the amount 
and velocity structure of the outer, unburned material. 

Some of the unburned material may be subject to a de- 
gree of downward mixing toward the inner ejecta when 
burning beco mes turbulent and enters the distributed 
flame regime (|Popelll987t INiemeven 119981 : lAspden et all 
2010). As of a result, it is not certain whether regions of 
unburned material are strictly located in the outer layer 
or mixed within the rest of the ejecta, nor is the mass 
range of unburn ed material from th ese thermonuclear ex- 
plosions known (jBaron et all 120031) . Therefore, the rate 
of detection of unburned material and its phenomenolog- 
ical details is of great importance for constraining how 
much material remains unburned in the explosion mod- 
els. 

The most prominent oxygen line in the optical spectra 
of SNe la is O I A7774. Unfortunately, since oxygen is a 
product of carbon burning, oxygen absorption lines are 
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likely to be a biased reference for measuring the location 
of unburned material. Thus one must look to carbon 
absorption features as a tracer for unburned material. 
While it is no t uncommon for ox ygen to be present in 
SN la spectra ([Branch et al.ll2006l ). signatures of carbon 
have only been sporadically reported in SN la spectra, 
most often during the pre-maximum phase. 

For the typical ejecta temperatures seen in SNe la 
(~10,000 K), the dominant i onization state of carbon 
is C II (jTanaka et all I2008D . At this temperature, 
the strongest optica l line of C II is that of A6580 
(jHatano et al.1 fl999b| ). When seen, this line produces 
a blue-shifted absorption near 6300 A that sometimes 
blends with the neighboring emission component of the 
Si II A6355 P-Cygni profile, thus making C II identi- 
fications problematic. There are wel l observed SNe la 



where C II A6580 is clearly present dPatat et alJ 


1996; 


Mazzali 2001; iGaravini et al.l 120051: Ifficken et al.l 


2007; 


Thomas et al.ll2007: Yamanaka et al.l2009b; Scalzo et al. 


2010), however in general, C II absorption is 


often 



weak and/or bl ended and therefo r e not a consp i cuous 
SN la feature d Salvo et al.l [20011; IBranch et ail l200l 
IStanishev et al.ll2007t) . ~ 



A number of factors intrinsic to the explosion it- 
self may contribute to the strength of C II A6580 and 
whether or not carbon features show up in SN la spec- 
tra. These include the asymmetrical distribution of car- 
bon, the extent of carbon burning, the temperature of 
the carbon-rich region, and the possible formation of an 
envelope of un burned material though a merger of two 
white dwarfs (iThomas et al.l l2007t iTanaka et all 120081; 
iScalzo et alj|2010f) . Low signal-to- noise spectra and ef- 
fects of line formation may also obscure weak C II A6580 
absorption features, thereby complicating the rate of de- 
tection further. 

In an attempt to better understand spectroscopic car- 
bon signatures in the pre-maximum optical spectra of 
SNe la, we present a comparative study of 16 SNe la 
where a C II A6580 absorption signature is evident, plus 
three additional cases where it may be present. Then us- 
ing the spectrum synthesis code, SYNDW, we fit this sam- 
ple of SN la spectra thereby mapping the velocity distri- 
bution of carbon-rich regions that give rise to C II A6580 
absorption features. To ensure a consistent analysis, we 
generated synthetic spectra for each of the time-series 
spectra in our sample. 

The new and archival data used in this study are pre- 
sented in §2. We review SN la subtype classifications in 
§3, and in §4 we discuss the spectrum fitting methods 
implemented with the SYN0W model, followed by our re- 
sults in §5. We discuss our results in §6 in the context of 
recent findings on the diversity of SNe la and conclude 
in §7 with our results regarding of the nature of C II 
absorption features. 

2. DATA 

Low signal-to- noise (S/N) SN la spectra make inves- 
tigations of any suspected C II absorption features dif- 
ficult. Therefore, with the possibility of carbon features 
being both weak and blended with Si II A6355, well ob- 
served multi-epoch confirmations of any suspected C II 
A6580 signatures were preferred for analysis. 

In constructing our data sample, we have included only 
SNe la where at least two consecutive spectra showed the 



presence of a ^6300 A absorption signature, indicating a 
likely C II A6580 absorption feature. In addition, if the 
C II feature was not seen before maximum light, then 
we did not include those SNe in our sample since this 
region of the spectrum becomes contaminated by neigh- 
boring Fe II lines w ithin a week after maximum light 
([Branch et al.l [20081 ). The two exceptions are the sin- 
gle epoch spectra of PTFlOicb and the 2002cx-like, SN 
2008ha, where the C II A6580 signature is evident. 

2.1. New Observations 

Pre-maximum optical spectra of SN 2010Y, 2010ai, 
and PTFlOicb show absorptions likely due to C II A6580 
and thus we included these in our sample. These data 
were reduced using standard IRAF procedures, were cor- 
rected for host galaxy redshift and are presented in Fig- 
ure Q] and Table 1. We used the s pectrum-comparison 
tool, SNID ([Blondin fe Tonrvll2007l ). to estimate the age 
of each spectrum relative to maximum light. 

Three spectra for SN 2010Y covering days —7, —6, and 
—3 with respect to maximum light were obtained using a 
Boiler & Chivens CCD spectrograph (CCDS) on the 2.4 
m Hiltner telescope at the MDM Observatory on Kitt 
Peak, Arizona. These spectra show the emergence of a 
relatively weak yet persistent 6330 A absorption feature. 
Two additional spectra were taken on day —2 and +1 
and were acq uired with the 9.2 m Hobby-Eberly Tele- 
scope (HET) ([Ramsey et al.ll 19981 ) at the McDonald Ob- 
servatory usi ng the Marcari o Low-Resolution Spectro- 
graph (LRS, IHill et al.lll998l ). By day -2, the 6330 A 
feature peaked in intensity and began to fade by day +1. 

Low resolution optical spectra of SN 2010ai were also 
obtained with HET. Spectra were taken well before max- 
imum light on days —10 and —8. These spectra show 
signs of a flux depression in the emission component of 
the Si II A6355 P-Cygni profile, suggesting the presence 
of C II A6580 absorption. 

Additional data comes from the di scovery of PTFlOic b 
by The Palomar Transient Factory ([Nugent et alJ feOlO). 
A low-resolution follow-up spectrum was obtained with 
the LRS on HET on June 3. The spectrum-comparison 
tool, SNID, identifies the spectrum of PTFlOicb as that 
of a normal SN la near day -10. Similar to SN 2010Y 
and 2010ai, the spectrum of PTFlOicb also exhibited a 
6330 A absorption feature. 

2.2. Archival Data 

In Table 2, we list these three recent objects along 
side 65 SNe la found in the literature with pre-maximum 
or near maximum light spectra. These 68 SNe have 
been o rganized by the subtype scheme of iBenetti et al.1 
(12005D and are labe led by the classification subtypes of 
IBranch et al.1 ([20061 ) (see §3). Since the progenitor chan- 
nel and the origin of C II A6580 absorption features 
may differ for Super-Chandra SNe la, we separate SN 
2003fg, 2006gz, 2007if and 2009dc as possible Super- 
Chandra candidates. We have also grouped SN 2000cx, 
2002bj, 2002cx, 2007qd and 2008ha as miscellaneous ob- 
jects since their cl assifications as SN la-like events are 
still under debate (IValenti et alJ[2059t iMcClelland et al.1 
[2011 [Poinanski et al.ll2010D . 

Optical spectra for all but three SNe la in our sam- 
ple were obtained from the Supernova Spectrum Archive 
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(SuSpeclQ) and sources therein. All spectra have been 
corrected for host galaxy re dshift and normalize d accord- 
ing to the formula given bv lJefferv etaLl (|2007ft in order 
to remove the underlying continuum since in this study 
we are only concerned with the position of the C II A6580 
absorption minimum and not absolute flux. The inclu- 
sion of SN 2008ha is for comparison purposes only and 
not meant to contribute to the discussion of carbon-rich 
regions of normal SNe la. 

In Figure [21 we present a single pre-maximum spec- 
trum of the 19 SNe la in our sample that show evidence 
of a C II A6580 absorption signature. Under the assump- 
tion of local thermodynamic equilibrium at 10,000 K, 
the strongest four optical lines of C II are AA4267, 4745, 
6580, and 7234, with the 6580 A line being the strongest. 
Rest frame positions of these lines are represented as ver- 
tical lines in Figure [5] and the four grey bands at the top 
indicate Doppler velocities of 1000-20,000 km s _1 . This 
corresponds to a region spanning 440 A wide and blue- 
ward of 6580 A, well within the extent of typical Si II 
A6355 P-Cygni profiles. 

In Table 2 we also note the epoch of the earliest spec- 
trum taken and how likely a C II A6580 detection is for 
each SN la (see §5.2). If both C II A6580 and A7234 
absorption features are present or if the C II A6580 is 
an obvious notch, then we denote the detection as "Def- 
inite" . If a weaker notch is seen, and both absorption 
wings remain, then the detection is labeled as "Proba- 
ble" . If a -6300 A feature comes in the form of a "slump" 
on the emission component of Si II A6355, then we note 
the C II A6580 detection as "Possible" . Some of the pub- 
lished data have S/N ratios that prevent making a clear 
case for C II A6580 detection and we note these as "Un- 
certain" . If none of these criteria are met, then there is 
"No" detection. 

3. SNE IA: SUBTYPE CLASSES 

The observed diversity of SNe la has been subdivided 
based on photomet ric a nd spectroscopic pr operties by 
IBenetti et all (|2005D and lBranch et all (|2006D . Below we 
briefly describe these classification schemes which will be 
r eferred to in our anal ysis . 

IBenetti et all (|2005l) grouped 26 SNe la according to 
two photometric and three spectroscopic observables, 
namely (1) Amis(-B) — the decline in m agnitude of th e 
B-band 15 days after B-band maximum (Phillips 1993), 
(2) Mb — the peak B-band magnitude, (3) iisui — the 
Si II expansion velocity rate of decrease, (4) i>io(Si II) 
— the expansion velocity of Si II A6355 10 days af- 
ter maximum light, and (5) lZ(Si) max — the ratio of 
AA597 2, 6355 absorption depth measured at maximum 
light (jNugent et al.lfl995l iBongard et al.lf2008h . 

With a sample of 26 SNe la, they found that SNe la 
could be organized into three discrete subtypes that were 
functions of velocity gradient and luminosity, namely (1) 
high- velocity gradient (HVG), (2) low- velocity gradient 
(LVG), and (3) FAINT. The subtypes HVG and LVG 
were mainly distinguished by having v values above or 
below —70 km s _1 day -1 , respectively, while the SNe 
in their sample with Mb max ~ —18.20 were labeled as 
FAINT. 

1 http://nhn.ou.edu/~suspcct 



iBranch et al.l (|2006ft took a purely spectroscopically 
bas ed approach classi fying 24 SNe la (later 65 SNe; 
see IBranch et al.l 12009) by the equivalent width of fea- 
tures near 5750 A and 6100 A which are usually at- 
tributed to Si II AA5972, 6355, respectively. When SNe 
la are arranged in this manner, the objects can be sub- 
divided roughly into four spectroscopic subtypes: (1) 
Core-Normal (CN) SNe la consist of objects such as 
SN 1994D, where from pre-maximum to 1-week post- 
maximum spectra are dominated by lines of Ca II, Fe II, 
Fe III, Mg II, O I, S II, Si II, and Si III, (2) Broad- 
Line (BL) SNe la are similar to CNs but instead display 
noticeably broader lines with higher average Doppler ve- 
locities, (3) Cool (CL) 1991bg-like spectra which exhibit 
low ionization energy ions such as Ti II, along with an 
increased ratio between the 5750 and 6100 A features, 
and (4) Shallow Silicon (SS) 1991T-likes display mostly 
high ionization energy ions, such as Fe III, in their pre- 
maximum spectra and are accompanied by weak S II 
absorption features. 

We note that comparing the sample of SNe la used in 
both studies shows that FAINT and HVG objects are 
equivalent to CL and BL respectively, and the LVG ob- 
jects are equivalent to CN and SS subtypes ([Branch et al.1 
I2006T ). Since we are using subtype classes that are mostly 
based on spectroscopic properties of SNe la, our study 
cannot address issues regarding broadband photometric 
observations. A systematic study of both the spectro- 
scopic and light-curve properties is difficult because these 
data can come from multiple sources and not every target 
we study here have both good spectroscopy and photom- 
etry. This type of s tudy is better suited for the Nearby 
Supernova Factory (lAldering et aJJ 12002ft and the Palo- 
mar Transient Factory (|Rau et al.ll2009l ). where the data 
sets are large and have good time-series coverage. 

4. SPECTRUM ANALYSIS MODEL: SYN0W 

One way to infer the velocity range of unburned ma- 
terial is to measure the absorption minimum of the C II 
A6580 line. The observed minimum of this feature is of- 
ten located near the strong emission component of Si II 
A6355, and is therefore subject to line blending and any 
limb- brightening that the Si II lin e contributes to the 
integrated spectrum (|Hoflich|[l990l ). Consequently, the 
observed minimum may underestimate the actual expan- 
sion velocity of the carbon-rich region. Therefore, in or- 
der to accurately estimate the true minimum of a blended 
line profile, one must take into account the effects of line 
formation by reconstructing the spectrum via numerical 
calculation. 

While there have been several new and more de- 
tailed spectrum synthesis codes presented in the liter- 
ature since the inception of SYN0W nearly thirty years 
ago, the SYN0W spectrum synthesis model remains a use- 
ful tool for the quick analy s is of resonant scattering line 
profiles dBaron et al.l Il994t iMazzalil 120001: [Branch et al.l 
120071: iKasen et al.l 120081: iBaron et all 120091) . Therefore, 
we chose the less computationally intensive approach of 
SYN0W to produce results that were internally consistent 
when we compared the SNe la in our sample. 

From a single-epoch optical spectrum, one can repro- 
duce many of the conspicuous features seen in SNe of 
all subtypes using SYN0W. However, to have greater con- 
fidence in the identification of a spectroscopic feature, 
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it is best to have a time series of closely spaced spectra 
in order to follow the evolution of both the observations 
and the fit. Because we wish to probe the nature of C II 
A6580 features in SNe la, we required fitting all spectra 
for the duration of time in which the 6300 A features 
could be seen. 

In SYNOW, one computes a spectrum by specifying the 
location and optical depth for a given set of ions. This 
allows one to infer spectral line identifications by directly 
fitting to a series of observed spectra. While SYNOW docs 
not calculate relative abundances of various elements, it 
is instructive at reconstructing the complex spectroscopic 
profiles brought about by multiple line scattering that is 
inherent to moving media. 

The version of the SYNOW model that we used can be 
described as follows: (1) a spherically symmetric and 
homologously expanding ejecta is modeled using a v oc 
r law, (2) light is emitted from a sharp photosphere, 
(3) optical depth, r, is a function of velocity as either 

(v/v p hot)~ n , exp »« , or exp ^ where each 
of these functions are characterized by indices n, v e , and 
a respectively, (4) line formation is purely due to reso- 
nant sca ttering and is treated using the Sqbolev Approx- 
imation (|Sobolevlll957l Uefferv fc Branchlfl990T) . and (5) 
for a given ion a reference line profile is calculated for a 
given t and the remaining lines follow from Boltzmann 
statistics. Input parameters for a SYNOW spectral fit in- 
clude: (i) a photospheric velocity (~v p hot)i (h) reference 
line r and minimum/ maximum velocities for each ion, 
and (iii) excitation temperature, T exc , to determine LTE 
level populations with respect to a reference line. 

As was done in the SN la comparative study of 
iBranch et afl (pOOl . for our invesitigation we too left the 
excitation temperature at 10,000 K for each ion. We also 
chose to use the exponential form of the optical depth 
profile with v e = 1000 km s . This aided in limiting 
the number of free parameters for a spectral fit. 

Because SYNOW assumes spherical symmetry, we are 
limited when investigating the possible asymmetrical dis- 
tributions of unburned material. If the distribution is 
that of a spherical layer, as is the case for the W7 model 
(|Nomoto et al.lll984t IThielemann et al.lll986D . then it is 
straightforward to make SYNOW model fits to compare 
with observations. In this case, the unburned material 
resides above the burned material in a spherical shell 
where the thickness depends on the extent of the burn- 
ing. For W7, this boundary is roughly at 14,000 km s _1 
with a stratified composition of IMEs and IPEs below. 

Recent multidimensional models of delayed detona- 
tions suggest that the unburned materi al may be left 
behin d in clumps throughout the ejecta (jGamezo et al.l 
120041) . For the situation where the unburned material is 
heavily concentrated to a single clump structure, we can 
only utilize SYNOW in certain cases (see §5.4). 

Observations have shown that there often exists higher 
velocity regions of line formation; namely Ca II and Si II 
(iHatano et all Il999at iKasen et all 120031 : iTanaka et all 
120101) . Such regions are said to be detached from the 
photosphere. We note a noticeable facet that arises when 
fitting Si II A6355 profiles during the earliest epochs is 
that the absorption width, the slow rise of the blue wing, 
and the sharp rise of the red wing require two separate 
velocity components of Si II to achieve a good match to 



observations. At best, a single-component of Si II with 
an increased value of the optical depth profile indices can 
only properly fit the red wing of the 6355 A absorption 
feature. 

G. H. Marion et al. (2011, in preparation) discuss re- 
cent observations of the Type la event, SN 2009ig, where 
there is clear evidence for both a photospheric region 
and a high velocity region of Si II. In their SYNOW analy- 
sis, they used the two-component approach in following 
the evolution of the Si II A6355 feature, which produced a 
better fit overall. Similarly, we too have adopted a proce- 
dure of using two components of Si II that are separated 
by ^4000—6000 km s _1 when necessary. 

One consequence that detached ions have on the line 
profile is that the emission component of the line is flat- 
topped. Some of our spectroscopic fits in §5.3 detach Si II 
to better fit the absorption component, while forfeiting a 
comparable fit to the full emission component; e.g. our fit 
for SN 1999ac. The impact that detaching has towards 
blending with C II A6580 is minimal and only one of 
offsetting the prescribed value of r for C II in the fit. 

5. RESULTS 

In terms of the spectroscopic diversity of SNe la, many 
of the differences are seen immediately after the explosion 
as the photosphere maps out the distribution of the out- 
ermost ejected material, with an i ncrease toward spec - 
troscopic conformity at later epochs ()Branch et all2 008) . 
Unfortunately, as shown in Figure El SN la spectroscopic 
observations during the two weeks prior to maximum 
light are underrepresented compared to those taken near 
maximum light or at later times. In addition, absorp- 
tion features attributed to C II A6580 are generally weak. 
As a result, they can be easily missed and thus go un- 
reported or are not securely verifiable due to low S/N 
and/or line blending. 

5.1. Frequency of C II Absorption Features 

Our search of the literature revealed that ^30% of SNe 
la with moderate to high S/N, pre- maximum spectra 
taken since 1983 January 1 show a feature near 6300 A 
that may be associated with C II A6580 absorption (see 
Table 2; Definite + Probable). Because of the sparsity 
of optical spectra before peak brightness, this percent- 
age may not represent the actual fraction of SNe la that 
exhibit C II A6580 signatures during the pre-maximum 
phase. Furthermore, the amount of published SN la data 
does not equal the amount of data taken, thereby placing 
considerable uncertainty on such an estimate. 

For comparison, in Figure [3l we also tally and plot 
the total number of SNe la where C II A6580 was first 
detected with respect to maximum light. When the to- 
tal number of published C II observations over the past 
couple decades is compared to the total number of SNe 
la discovered during the pre-maximum phase in a span 
of three years, the occurrence of absorption due to C II 
A6580 appears infrequent. 

However, many SNe la might show C II (or C I, or 
C III; see §6.1) during the first weeks of the explosion, 
but may not exhibit strong, identifiable carbon absorp- 
tion features by the time they are discovered and spectra 
taken. If spectra are not obtained prior to maximum 
light, this can lead to an under-reporting of C II SNe la 
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cases. For example, in some objects, the C II features 
appear to fade by day —5 or even earlier. This raises 
the question: Might all SNe la show appreciable C II ab- 
sorption at some level if observed at a sufficiently early 
enough epoch? 

Given the weak impact that C II features are observed 
to have on the integrated optical spectrum, identifying 
every instance for when C II is present rests heavily on 
one's ability to distinguish the signature from noise and 
the effects of line blending. In Figure |H we compare the 
spectrum of SN 1998aq at day —9 to three cases where 
we have artificially added Gaussian noise to the original 
spectrum. With only 10% and 20% noise added, the C II 
A6580 feature remains distinguishable. However, at 25% 
the feature begins to lose its identity above the noise. An 
example of this is t he spectrum of SN 2 05bl a t day —6, 
shown in Figure |H iTaubenberger et all (2008) proposed 
a C II A6580 identification for the 6400 A depression seen 
in the spectrum. However, we excluded this SN from our 
sample since we were not able to generate an accurate fit 
for the feature. 

5.2. SNe la C II Features: Conspicuous to Weak 

Like most features in the spectra of SNe la, the de- 
gree of adjacent line blending that C II A6580 absorption 
signatures undergo can obscure the full extent of a line 
profile. More often than not, the minimum of the C II 
A6580 line is blended with the P-Cygni emission compo- 
nent of the Si II A6355 line. The day —5 spectrum of SN 
1999by is a good example of when this takes place (see 
Figure EJ. 

When the C II A6580 absorption feature is weak there 
is also often no obvious C II A7234 signature. However, 
when the C II A6580 feature is strong, a corresponding 
C II A7234 feature does begin to appear as the excited 
level of the 7234 A line becomes more populated. There 
are several cases, such as SN 2006gz and 2007if, where 
C II A6580 is certainly present but either the wavelength 
coverage is incomplete or the spectrum is too noisy to 
say whether or not C II A7234 is present as well. 

In Figure El we plot two spectral regions of four SNe 
la in Figure [H that highlight cases when both of these 
C II lines are present. For each SN, the red and blue 
lines denote AA6580, 7234 Doppler velocity-scaled spec- 
tra, respectively. The absorption minima overlap nicely 
(see dashed vertical lines) and the symmetry about the 
minima is indicative of a match between the spectral sig- 
natures of the same ion. We interpret this to mean that 
the 7234 A line is indeed present when the 6580 A line 
is strong. Currently, these two C II lines are the best 
means by which to securely identify whether carbon is 
present in the early-epoch SN la spectra. 

The C II lines that appear further in the blue have 
also b een suggested to be present as well (jThomas et al.l 
2007). However, this region is crowded by lines of 
IPEs making it more difficult to securely identify C II 
AA4267, 4745 and use for determining accurate Doppler 
velocities. This appears to even be the case for SN 
2006gz, 2007if, a nd 2009dc where the C II A6580 ab- 
sorption is strong (IHicken et al J 120071 ; IScalzo et al J 12010 : 
ITaubenberger et al.ll2011ft . 

5.3. SYNOW Model Fitting 



For each of the 19 SNe la in our C II sample, we 
produced a time series of synthetic optical spectra that 
covered the observed extent of the C II A6580 feature's 
presence as well as the observed wavelength coverage of 
the data (~3500-9000 A). Our initial fits included the 
canonical set of IMEs and IPEs that are prevalent in 
pre-maximum SN la spectra, whereby afterward we in- 
cluded C II and adjusted the optical depth and detach- 
ment velocity until a match to observation was made. 
The goodness-of-fit to the observed Si II A6355 — C II 
A6580 blended profiles did not change with the full set of 
IMEs and IPEs removed. Thus, for each SN in Table 3, 
we only list the relevant parameters for C II and Si II. 

As was discussed above in §4, some of the observed 
Si II A6355 profiles required two separate components of 
Si II in order to fit the full width of the absorption. We 
indicate which SNe la in our sample required this fitting 
procedure by appending the inferred Doppler velocities 
and optical depths with ">" in Table 3. 

In Figure [6j we compare fits to a single observed spec- 
trum for each of the SNe la in our sample. The black 
lines represent the observed spectra, the red lines are 
SYNOW fits where C II has been included, and the blue 
lines are the same fits without C II. The synthetic spec- 
tra match fairly well to a variety of Si II A6355 — C II 
A6580 blended profiles, and the interpretation that the 
6300 A feature is due to C II A6580 is in agreement with 
that of previous authors (see references in Table 2). 

In this study, we also offer new and revised expansion 
velocity estimates of C II for a couple SNe la , particularly 
SN 19 90N and 1999by. It was suggested bv lFisher efall 
(|1997f) that the overly broadened 6040 A absorption in 
the day —14 spectrum of SN 1990N was due to a two- 
component blend composed of Si II at ~20,000 km s -1 
and C II at -26,000 km s" 1 . Similarly, iMazzalil (|2001f > 
suggested that this feature is predominantly due to Si II 
A6355 but also requires an outer zone of high velocity 
carbon between 19,000 and 30,000 km s" 1 . If this in- 
terpretation is correct, then other SNe la like SN 1990N 
would also require a similar zone of carbon to reproduce 
their early-epoch spectra. However, our fit for SN 1990N 
uses two components of Si II to fill the 6040 A feature 
while C II is only at 16,000 km s _1 to account for the 
6300 A feature. 

In the case of the s ub-luminous SN 1999by , C II A6580 
was not reported by IGarnavich et all (|2004l) . However, 
our fit in Figure [6] is fairly convincing when C II is in- 
cluded as a detached layer 2000 km s above the photo- 
sphere for the day —5 spectrum. In fact, almost 80% of 
the SNe in our sample suggest at least a mildly detached 
layer of C II during some point along the evolution of the 
absorption feature. For example, we modeled the C II in 
SN 1994D initially at 14,000 km s" 1 coincident with v p hot 
at day —11, after which it remains at this velocity as the 
photosphere recedes to 12,000 km s^ 1 by day —5. 

6. DISCUSSION 

The strength and velocity range of carbon in pre- 
maximum SNe la spectra can provide a valuable tool 
to investigate various explosion models. The W7 defla- 
gration model, for instance, contains a ~0.07M Q layer of 
unburned material above 14,000 km s _1 . Using the spec- 
trum synthesis and model atmosphere code, PHOENIX, 
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iLentz et"al"1 (|200lD compared calculated non-LTE spec- 
tra to spectroscopic observations of the normal la event, 
SN 1994D. Despite the outer layer of unburned mate- 
rial in W7, none of their synthetic spectra were able to 
account for the observed C II A6580 absorption feature 
near 6290 A. Our SYNOW fits for this object allow for C I I 
below 14,000 km s _1 . To that end. iTanaka et al.1 (|2011f ) 
were able to reproduce a C II A6580 absorption feature 
seen in the day —11 spectrum of SN 2003du by placing 
the carbon-rich region at lower velocities in their version 
of W7. They obtained an upper limit on the abundance 
of carbon to be 0.016-M© at v > 10,500 km s _1 in this 
object. 

At late-times and based on the t hree-dimensional de- 
flagration models of lRopkd (|2005[ ). IKozma et all (|2005Q 
compared synthetic spectra to the late-time spectra of 
three SNe la and set an upper mass limit of unburned 
material below 10,000 km s" 1 to be -0.07 M . Other 
estimates for the mass of unbur ned material have been 
made using delayed-detonations (Ho flich et al.H2002h and 
other modeling; e.g. lower limit of O.O14M of unburned 
material between 10 00 and 14,000 km s" 1 for SN 2006D 
(|Thomas et al.l |2007[ ). However, since each estimate is 
obtained by different means, a comparison of such re- 
sults does not advance the discussion on the nature of 
C II A6580 absorption features. 

In order to utilize C II A6580 absorption features for 
estimating the mass of unburned material, the effects of 
temperature and geometry of the carbon-rich regions as 
well as the influence of radiative transfer effects must be 
explored. Our SYNOW modeling begins this process by 
mapping the observed velocity distribution for a sample 
of 19 SNe la. Below, we discuss our interpretation of the 
observed frequency of C II absorption features and how 
they relate to the properties of carbon-rich regions and 
SN la diversity. 

6.1. Temperature Effects on Carbon Features 

The temperature of the carbon-rich region will influ- 
ence whether or not C II is the dominant ionization 
species and will therefore dictate the strength of C II 
A6580 a bsorption features. Using non-LTE calculated 
spectra, iNugent et al.l (1995) pointed out that the spec- 
troscopic sequence observed among various SNe la could 
be explained by a continuous change in the effective 
temperature of the ejecta (7400—11,000 K), from cool 
1991bg-likes to the hotter 1991T-likes (see their Fig. 1). 
While the ejecta of SNe la are an environment with 
non-LTE processes, under the assumption of a C+O-rich 
composition the Saha-Boltzmann equation indicates that 
carbon will mostly be in the form of C II between 6000 
and 12,000 K. 

Of the 68 SNe la listed in Table 2, 10 out of the 14 
CNs and five out of the 13 CLs exhibit C II A6580 ab- 
sorption features. From a nearly equal sampling of SN la 
subtypes, the fact that two CNs for every one CL SN la 
exhibit C II absorption lines suggests that the presence 
of C II absorption features depend on the effective tem- 
perature to some degree. This point is also consistent 
with the fact that only two of 13 SSs show signs of C II 
A6580 absorption features in their spectra. 

We note that C I absor ption features are no t often de- 
tected in SN la spectra. iMarion et a l. (2006) presented 
NIR spectra of three normal SNe la and discussed the 



lack of C I absorption signatures due to the absence of 
the strongest C I NIR lines, namely C I AA9093, 10691. 
A more extensive study of 41 SN la NIR spectra span- 
ning two weeks before a n d afte r maximum light was dis- 
cussed by IMarion et al.1 (|2009f ) and they too reported a 
lack of C I signatures. Because spectral signatures of car- 
bon burning products were observe d to occupy the same 
region of ejecta (Mg II and O I), IMarion et all (|2006f) 
concluded that nuclear burning had been complete out 
to at least 18,000 km s" 1 in their objects. 

It is perhaps not surprising that many NIR spectra of 
SNe la do not show lines from C I, given that (1) the 
abundance of carbon may be only ~1-10% of the total 
ejected mass, and (2) if there is carbon present in the 
outer lay ers then it is mostly once ionized. Interestingly, 
however. iHoflich ~et al. (2002) reported a conspicuous C I 
A10691 absorption feature in the day —4 NIR spectrum of 
SN 1999by. They were able to reproduce many features 
of the observed optical and NIR spectra using a series of 
sub-luminous delayed-detonation models with a range of 
transition densities, pt r , between 8—27 x 10 6 g cm~ 3 . 

The appearance of the C I A10691 line was seen con- 
currently with the optical C II A6580 absorption feature 
we reported above in §5.3. Both the C I and C II spec- 
tral signatures indicated that the carbon-rich region was 
above the 12,000 km s _1 photosphere; i.e. the minimum 
of the C I A10691 absorption corresponds to —13,000 km 
s _1 while our SYNOW fits for this object place C II at 
14,500 km s — At least for this cool sub- luminous SNe 
la, the influence of temperature on the ionization state of 
the carbon-rich region is apparent from the simultaneous 
appearance of C I and C II absorption features. 

Additional evidence that ejecta temperature plays a 
role in the detection of carbon in SN la spectra would be 
if C III were clearly detected in a hotter SN la subtype, 
such as SN 1991T or SN 1997br. The similarity between 
the ionization potentials of C II (24.4 eV) and S II (23.3 
eV) suggests that the presence of C III absorption fea- 
tures may be concurrent with spectroscopic signatures 
of S III. Fortunately, we can examine optical spectra to 
check for the simultaneous presence of C III A4649 and 
S III A4254 absorption features. 

In Figure we plot and compare SYNOW fits for 
SN 1991T and 1997br where we have included C III 
and S III. The identification of the 4500 A feature of 
SN 19 91T has been discussed before b y lHatano et al. l 
(2002) , and similarly in other SN la (jGaravini et al.l 
20041: iChornock et all 120061 ). While it was argued by 
Hatano et al.l (|2002l ) that including C III in the fit pro- 
duced a mis-match with the observed spectrum (too blue 
overall) , our synthetic spectra (red lines) are in fair agree- 
ment with observations near 4500 A for both SN 1991T 
and 1997br. 

In addition, a better fit to the 4250 A absorption fea- 
ture is obtained with the inclusion of S III. The 4250 A 
feature is predominately due to the Fe III A4404 multi- 
plet. However, by adding S III to the fit we were able to 
fill in the blue wing of this absorption feature for both 
objects. Therefore, our identification for the 4500 A ab- 
sorption feature as being that of C III A4649 is more 
likely, though the evidence is circumstantial. If the C III 
A4649 identification is correct, then this may indicate 
that a lack of C II (or C I) spectroscopic features do not 
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necessarily imply the complete burning of carbon in the 
hotter subtype events. 

6.2. Interpreting C II X6580 Doppler Velocities 

In this paper, our attention has been primarily focused 
on carbon-rich regions in the outermost ejecta. The wide 
range of observed C II Doppler velocities among different 
SNe la suggests a large variation in the extent of carbon 
burning, with some objects exhibiting high velocity car- 
bon while in others the carbon is present at low velocities 
(see Table 3). How high or low is usually in reference to 
the position of the carbon cut-off seen in the W7 model 
(~14,000 km s _1 ) instead of the kinetic energy of the 
supernova itself. 

6.2.1. Extent of Burning via C II \6580 

Given that the characteristic ejecta velocity is propor- 
tional to (Ekin/Mej) 1 / 2 (|ArnettlH98l ). a standardized 
way of looking at the extent of burning is a more appro- 
priate measure for interpreting the range of observed car- 
bon velocities. The 6355 A line of Si II has been used as 
an indicator of the photospheric velocity at early epochs 
(jJefferv fe Branch! fl990t iPatat et al.l fl996h . Choosing 
Si II A6355 over other lines alleviates any difficulty in 
obtaining Doppler velocities amid too much line blend- 
ing and allows for consistent time-coverage. This makes 
the absorption minimum of Si II A6355 a good point of 
reference for investigating the extent of burning via C II 
A6580 absorption features before maximum light. 

In Figure we have the ratio of Doppler velocities, 
v(C II A6580)/z;(Si II A6355), derived from our SYNOW 
fits plotted versus days relative to maximum light. This 
shows that (1) for an individual SN, the ratio remains 
at a fairly sustained value over time and (2) the differ- 
ent velocity ratios among the SNe lie roughly within the 
same region and are similar to within ±10%. We have 
ignored the three outliers because one is in a region of 
the plot where line blending obscures the supposed C II 
A6580 detection (SN 2001V; see §6.2.2) and the other 
two objects suggest that the carbon is clumpy and not 
along the line-of-sight of the observer (SN 2006bt and SN 
2008ha). 

The notion of an optically thick photosphere in SNe is 
generally a good assumption, even though the line form- 
ing region may extend 500 km s _1 in either direction. 
Therefore, for a given spectrum, any velocities of an ion 
that are measured to be fairly below that of the pho- 
tospheric velocity, v p hot, may indicate ejecta asymme- 
tries. That is, any observed discrepancy between v p hot 
and veil could be explained if the actual velocity of C II 
is the same as v p hot but instead of forming in a shell at 
the observed velocity, the carbon is in a clump at v p hot 
and offset by an angle, 0, from the line of sight. In par- 
ticular, we can estimate this projection angle if van < 
v phot . For SN 2006bt and SN 2008ha we calculate this 
projection angle to be ~50° and 60°, respect ively, where 
our resu lt for SN 2006bt is in agreement with lFolev et al.l 

(poiobl) . 

If the dominant C II behavior is due to asymmetrical 
distributions of a single localized clump of unburned ma- 
terial, one could expect there to be more scatter below 
the mean that is presented in Figure [8] since an arbitrary 
orientation of the clump relative to the observer ought to 



lead to more cases where i>(C II A6580)/t>(Si II A6355) < 
1 (like SN 2006bt). Instead, what we find for the SNe in 
our sample is that the distribution of carbon-rich mate- 
rial is consistent with a layered or hemispheric geometry. 

Another surprising aspect of Figure |8] is that the can- 
didate Super-Chandra SNe la resid e in the same region 
of the plot as the other objects. iScalzo et al.l ()2010f ) 
suggested that the large C II A6580 feature, concurrent 
with low velocities, could be explained by invoking a 
pre-explosion envelope of progenitor material originat- 
ing from the merger of two white dwarfs. In this sce- 
nario, the explosion is inhibited by and loses kinetic en- 
ergy to the envelope, ionizing the shell of surrounding 
carbon. Whatever the nature of carbon-rich regions in 
these three SNe la, it is also constrained by the value of 
v(C II A6580)/w(SiII A6355). 

6.2.2. Consequences of Line Blending 

Toward the uppermost region of Figure [H there is a 
noticeable lack of highly detached C II objects with time- 
series coverage that exhibit pre-maximum C II A6580 ab- 
sorption features. This does not necessarily imply that 
carbon-rich regions in SNe la are absent above a par- 
ticular velocity. Rather, the "missing" subset of objects 
may be a result of SNe la with lower density carbon-rich 
regions further out and/or radiative transfer selection ef- 
fects; e.g. line blending. 

Most of the SNe la in our sample that exhibit C II 
A6580 absorption features are of the LVG subtype while 
only one is an HVG event (SN 2009ig). This is either a 
real trend of C II A6580 absorption features and there- 
fore a diagnostic of SN la diversity, or a result of line 
blending due to very high velocities of carbon-rich re- 
gions. In regards to the latter cause, some SNe la might 
contain a C II A6580 absorption feature in their spectra 
but its presence is obscured by the Si II A6355 absorption 
trough, particularly for HVGs. 

Assuming that the scatter of v(C II A6580)/v(Si II 
A6355) values in Figure [5] is the same for both HVG and 
LVG subtypes, and assuming that HVGs span a v p hot- 
space from 14,000 to 16,000 km s _1 , then the minimum 
of any C II A6580 signature would be between 6120—6280 
A. If weak and blue-shifted to these wavelengths, the C II 
A6580 absorption would most likely go undetected. If the 
C II optical depth were large enough, then the most di- 
rect evidence for hidden C II A6580 would be whether or 
not an associated absorption from C II A7234 could be 
seen between 6730-6900 A. 

This raises the question: At what velocity will C II 
A6580 completely blend with Si II A6580? In FigureEJ we 
presented observed C II A6580 blending scenarios where 
the signature is weak and nearly hidden, such as possibly 
observed in SN 2001V and 2009ig. These objects suggest 
that C II A6580 absorption features can in some cases 
be obscured via line blending. Such a "hidden" C II 
signature would be seen in the form of an absorption on 
the shoulder on the blue wing of the Si II A6355 emission 
component. 

A series of synthetic spectra that include only blends 
of Si II A6355 and C II A6580 line profiles are shown in 
Figured] (top panel). For these SYNOW spectra, the veloc- 
ity of Si II was fixed at 10,000 km s" 1 while the velocity 
of C II was increased from this 10,000 km s _1 by 1000 
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km s -1 increments up to a velocity of 15,000 km s -1 . 
We simultaneously decreased the optical depth in keep- 
ing consistency with the observed profile shapes. At a 
velocity of 15,000 km s _1 , C II A6580 loses its discerni- 
bility as a feature and remains hidden until C II reaches 
a velocity of 27,000 km s _1 . At present there are no 
observations reporting a C II A6580 feature appearing 
blue-ward of the Si II A6355 absorption. 

Aside from being able to account for the observed va- 
riety of C II absorption profile shapes and velocities with 
respect to Si II A6355, in Figure [9] we also show that the 
shoulder effect can be reproduced when C II is detached 
with a velocity —4000 km s _1 greater than that of Si II. 
This is consistent with our fit for SN 2001V where we 
have C II and Si II at 12,000 and 8000 km s _1 , respec- 
tively. 

In Figure [9] (bottom panel) , we show a series of syn- 
thetic spectra where C II and Si II are placed at the same 
velocity and increased in sequence from 10,000 to 18,000 
km s _1 . We modeled the Si II profiles using the two- 
component method so as to reproduce HVG-like proper- 
ties. As can be seen, the C II A6580 absorption feature 
can be easily obscured at higher velocities. The proposed 
weak C II A6580 feature for SN 2009ig lies atop the emis- 
sion component of Si II A6355 and is consistent with this 
series of synthetic spectra. 

It is generally thought that the infrequent number of 
HVGs where C II A6580 is detected suggests an environ - 
ment of sparse carbon-rich regions (jPignata et al.ll2008h . 
However, if the proposed C II A6580 signature in SN 
2009ig is correct, then our fits suggest that the feature 
could be weak based on line blending, and not necessar- 
ily due to a low carbon abundance (or lower than that of 
LVGs) alone. 

6.3. Asymmetries in the Distribution of C II Regions 

Because the ejected material follow a homologous ex- 
pansion law (v oc r) , the degree to which electron scatter- 
ing will polarize emergent light across a spectral line can 
provide geometrical information o n the material between 
the o b server and the photosph ere (|Shapiro fc Sutherland! 
119821: IWang fc Wheelerl 12008ft . Spectropolarimetry of 
some SNe la have revealed that the distribution of 
the outermost ejected material has an ov erall deviation 
from spherical sym metry of up to —10% (Hoflich| ll991l ; 
IWang et al.l 11997ft . Specifically, optical Si II P-Cygni 
profiles are observed to have peak polarizat ion values of 
0.3— 2.0%, five days before maximum light (Wa ng~e"t al.l 
2007). Other absorption features, such as the Ca II IR 
triplet, have also been observed to be strongly polarized 
(Wa ng et al.ll200l IKasen et al.ll2003l) . 

Of the 19 SNe la in our sample, only SN 1994D, 
1999by and 2005hk have nearly simultaneous spectro- 
scopic and polarization data during the presence of a 
C II A6580 absorption feature. Polarization levels for 
these objects across Si II A635 5 were reported as insignif- 
icant (-0.3% for SN 1994D. IWang et al.l [19961: -0.4% 
for SN 1999bv. IHowell et al.l l2001t -0.4% for SN 2005hk, 
iChornock et al.ll2006l ). whereas the degree of polarization 
was even less in the region of the C II A6580 signature. 

Regarding S N 1999by, th e thre e spectropolarimetric 
observations of iHowell et al.l (|2001l ) were combined from 
data taken on days —2, —1 and in order to increase 
the S/N ratio (see their Fig. 2). The resultant spectrum 



exhibited the same 6300 A absorption feature that we 
attrib uted to C II A6580 in the day —5, —4 and —3 spec- 
tra of iGarnavich et al.l (120041) . but t he concurrent po- 
larimetrv data of IHowell et al.1 (|2001ft did not show any 
significant amount of polarization near the feature. This 
would suggest that the C II regions are roughly spherical 
in at least these objects. However, because spectropo- 
larimetric observations are most useful when a spectral 
line is strong and the S/N is high, it is difficult to use 
the weak C II A6580 absorption features of SNe la to in- 
fer asymmetries of carbon-rich regions with polarization 
data. Despite this, possible asymmetries of carbon-rich 
regions mig ht be gleaned throug h other means. 

Recently, IMaeda et all (|2010bh discussed SN la diver- 
sity in the context of global asymmetries of ejected ma- 
terial. They argued that the HVG and LVG subtypes 
constitute a picture of SNe la in terms of an off-center 
delayed-detonation at different viewing angles. In this 
scenario, an off-center ignition is followed by the propa- 
gation of a sub-sonic deflagration flame that imprints an 
off-set distribution of high density ash. Once the transi- 
tion to detonation occurs, a super-sonic flame only suc- 
cessfully burns material in the lower density regions and 
is effectively screened by some of the deflagration ash. 
What remains is a lopsided distribution of burning prod- 
ucts, resu lting in a hem i spheri c asymmetry of the ejecta. 
Similarly, Ma und et al.l (|2010l ) reached the same conclu- 
sion regarding SN la diversity after examining a possible 
relationship between HVG /LVG subtypes and the possi- 
ble asymmetrical distribution of photospheric Si II. 

In this model, LVGs correspond to viewing the hemi- 
sphere of ejecta that coincides with the off-center igni- 
tion, whereas HVGs are the contrary with an opening 
angle of —70—75°. Since we do not see a predominate 
number of HVGs with conspicuous C II A6580 features 
in our sample, this is consistent with the above model, 
owing to the fact that detonation waves ought to leave 
little unburned mat erial behind on the HVG side. If 
IMaeda et al.l (|2010bft is correct, and if our sample of SN 
la that exhibit C II A6580 is representative, then this sug- 
gests that the filling factor of carbon in the LVG hemi- 
sphere is less than unity. 

7. CONCLUSIONS 

In an effort to better understand the frequency and 
general properties of C II absorption features, we exam- 
ined C II A6580 signatures in the pre-maximum spectra 
of 19 SNe la, of which included 14 "normal" SNe la, 
three possible Super-Chandra SNe la, and two 2002cx- 
like events. Using SYN0W to produce synthetic spectra, 
we modeled observed —6300 A absorption features as a 
C II A6580 P-Cygni profile blended with that of Si II 
A6355. Through our SYN0W model fits we estimated the 
C II expansion velocities for a variety of objects. Below 
is a summary of our major findings: 

1. A survey of the optical spectra of 68 objects pub- 
lished in the literature since 1983 January 1 indi- 
cates that up to —30% of SNe la display an absorp- 
tion feature near 6300 A that can be attributed to 
C II A6580. While this percentage is likely biased, 
it does suggest that C II A6580 absorption features 
are more common t han was previously suspected 
(jThomas et al .1120071) . If spectroscopic observations 
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of SNe la are obtained more than ~1 week before 
maximum light, we suspect an even larger fraction 
of SNe la of all subtypes may show identifiable C II 
absorption signatures. 

2. A greater frequency of C II A6580 absorption fea- 
tures appear in the low-velocity gradient subtypes 
(LVGs) compared to high-velocity gradient events 
(HVGs). Th i s is in line with the interpretation of 
iMaeda et al.l (|2010bf ) , supporting the idea that part 
of SN la diversity can be accounted for by viewing 
angle and off-center ignition effects. 

3. The influence of the temperature of the carbon-rich 
region on the incidence of C II A6580 signatures is 
most evident for the 14 "normal" SNe la in our 
survey. The frequency of C II A6580 signatures 
peaks with the Core-Normal (CN) SNe la, while 
the Cool (CL) and Shallow-Silicon (SS) subtypes 
exhibit fewer C II A6580 detections. This result is 
co nsistent with t he eff ective temperature sequence 
of lNueent etal] (fl995l) . 

4. We find the values of v(C II A6580)/^(Si II A6355) 
among 16 of the SNe la in our C II sample are sim- 
ilar to within ±10%. Assuming the minima of Si II 
A6355 absorption features are an appropriate mea- 
sure of photospheric velocities prior to maximum 
light, then the small number of cases where v(C II 
A6580)/t>(Si II A6355 < 1 could be indicative of ei- 
ther a layered distribution or multiple clumps with 
a comparable filling factor. 

One of the most interesting results of this study is that 
C II A6580 absorption features might be in many, if not 
most, early-epoch SN la spectra. We initially set out 
to investigate what we suspected was a small number 
of known detections of C II A6580. However, when we 
examined the 6300 A region in a relatively large sample of 
objects, we discovered that weak signatures of C II A6580 
were often detectable and could influence the observed 
Si II A6355 spectral profile. 

One way to test for a high frequency of C II A6580 fea- 
tures in pre-maximum SNe la spectra would be to obtain 
high S/N observations of the region between 6700—7200 
A in order to look for the relatively weaker C II A7234 
line that should also be present. In addition, a data set 
of C II A6580 observations with wavelength coverage that 
encompasses the O I A7774 absorption feature would al- 
low for investigating any correlations between C and O 
spectral line properties. 

Understanding the presence of C II in SN la spectra 
may prove to be a valuable diagnostic of SN la diversity 
and therefore a possible means by which to probe the 
underlying explosion mechanism. Utilizing the absolute 
strength of the C II A6580 line in order to extract abun- 
dance information will require comparison to synthetic 
spectra calculations based on hydrodynamical modeling. 
A detailed comparative study that can extract C abun- 
dance information, derived from C II A6580 absorption 
features, promises to produce a wealth of insight regard- 
ing carbon-rich regions, and therefore constrain the pa- 
rameter space of hydrodynamical modeling. 
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Figure 1. Pre-maximum spectra of SN 2010Y, 2010ai, and PTFIOicb. Data are from HET and the 2.4 m Hiltner at MDM, have been 
scaled appropriately, and are in the rest frame of their respective host galaxy. Absorption features that are consistent with C II A6580 
blue-shifted to typical expansion velocities of SNe la are seen at ~6300 A. See Table 2 for observational details. 
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Figure 2. Sample of 19 SNe la whose pre- maximum optical spectra contain an absorption due to C II A6580. A pre-maximum spectrum 
of SN 2004dt (HVG) is shown at the bottom (dashed line) for comparison. Spectra have been normalized and corrected for host galaxy 
redshift. Vertical lines indicate the rest wavelength position of C II AA4267, 4745, 6580, and 7234. The width of the shaded regions near 
each vertical line corresponds to Doppler velocities between 1000 - 20,000 km s — . To guide the eye, a couple black arrows indicate the 
position of a likely absorption feature of C II A6580. 
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Figure 3. Shaded in grey are the number of SNe la discovered between 2006 January 1 and 2009 December 1 plotted against the epoch 
at which they were identified as being a Type la supernovae. Shaded in black are the number of SNe la that show ~6300 absorption 
features in their pre-maximum spectra at the time when the feature was first seen. Data were taken from Central Bureau Electronic 
Telegrams (CBET) and Astronomical Telegrams (AT EL), Estimates for the age of each sup ernova were o b tained by one of several publicly 
availa ble spectru m-comparison tools ( PASS parToo; [Harutvunvan ct al. (2005), Superfit; IHowell et all (2003), SNID; IBlondin fc Tonrvl 
lf2007ft . GELATO; Harutv unvan et al.l lf200l )1. The cited epoch for each SNe la in this plot was taken as reported in the telegrams. The 
numbers above are not a strict sampling of all objects during the four year period, given that some of the age estimates were too vague to 
be included. This explains the peaks at day — 7, 0, +7 and +14 relative to maximum light. Because the spectral identification programs 
determine age with several days of error, we binned the data sample at days —13 ± 1, —10 ± 1, —7 ± 1, —4 ±1, + 2, +4 ± 1, +7 ± 1, 
+10±1 and +13 ±1, Both histograms show, on average, that observations before 1— week pre-maximum are not well represented compared 
to observations near maximum light or later. 
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Figure 4. Five segments of spectra showing the questionable of a C II A6580 absorption feature as a function of the signal-to-noise ratio. 
At the top is the original spectrum for the day —9 observation of SN 1998aq. The three spectra that follow represent the day —9 spectrum 
with 10, 20, and 25% Gaussian noise added. The spectrum at the bottom is that of SN 2005M, six days before maximum light, where 
ITaubenbereer et al.1 <205sft identified the 6400 A feature with that of C II A6580. 
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Figure 5. Doppler velocity scaled spectra for four of the SNe la in our sample. Blue and red lines represent C II A6580, 7234, respectively, 
for each SN. Dashed vertical lines are positioned at the minima of the absorption features to guide the eye. 
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Figure 6. SYNOW fits of the 19 SNe la in our sample presented in two panels. The red and blue lines denote a synthetic fit to the spectrum 
of a SNe la with and without C II, respectively. The contrast between both fits suggests the presence of C II with varying strength in these 
objects. See Table 2 for fitting parameters. 




Figure 7. SYNOW fits for SN 1991T and 1997br, two SS SNe la, around 1-week pre-maximum. In addition to matching the major features 
with Fe III, S II, Si II, and Si III, the absorption at 4500 A is identified with a multiplet of C III A4649. As before in Figure [6] the red line 
denotes the synthetic spectrum with carbon, and the blue line is without. 
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Figure 8. v(C II A6580)/i)(Si II A6355) vs. epoch for our C II— SNe sample. Circles show CNs, side- ways diamonds show BLs, squares 
show CLs, diamonds show SSs, upside down triangles show Super-Chandra candidates, and the 'X' denotes the 2002cx-like, SN 2008ha. 
Calculated error assumes spectral fits to the velocity of an ion are good to within 500 km s . The black horizontal line represents the 
average value of data points (excluding SN 200IV, 2006bt and 2008ha) and the dashed horizontal lines indicate a 10% difference from the 
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Figure 9. (top) A series of synthetic spectra show an evolution of blending scenarios between C II A6580 and Si II A6355 as the velocity of 
C II is increased from 10,000 to 15,000 km s _1 . The two spectra at the bottom represent a case where the C II is present at high velocities. 
Note that even if the 6580-line is obscured, the 7234-line is blueshifted to the rest frame position of the 6580-line. Three high-velocity 
scenarios are also shown to note the velocity at which C II A6580 emerges from the blue wing of Si II A6355. (bottom) Like the top figure, 
but instead C II and Si II are at the same velocity and increased in sequence together while holding t fixed. The Si II profiles are made 
from two-components of Si separated by 4000 km s — 1 . The day —14 spectrum of SN 2009ig is shown for comparison. 
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Table 1 

Summary of Observations 



Supernova 


Host Galaxy 


Rcdshift 


Observation 


Epoch 


Telescope 


Name 




(km s- 1 ) 


Date 


(days) 


Instrument 


SN 2010Y 


NGC 3392 


3256 


2010 Feb 09 


-7 


MDM/CCDS 








2010 Feb 10 


-6 


MDM/CCDS 








2010 Feb 13 


-3 


MDM/CCDS 








2010 Feb 14 


-2 


HET/LRS 








2010 Feb 17 


+1 


HET/LRS 


SN 2010ai 


SDSS J125925.04+275948.2 


5507 


2010 Mar 11 


-10 


HET/LRS 








2010 Mar 13 


-8 


HET/LRS 


PTFlOicb 






2010 June 3 


-10 


HET/LRS 
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Table 2 

SN la Sample 



SN 
Name 



Sub-type a 



C II? 



Epoch Spectrum 
(days) Source 



SN 
Name 



Sub- type 



C II? 



Epoch 
(days) 



Spectrum 
Source 



LVG 



SN 1990N b 


CN 


Definite 


-14 


M01 


SN 1999dq 


SS 


No 


-9.5 


M08 


SN 1991T b 


SS 


Possible 


-13 


F99 


SN 1999ee 


SS 


Uncertain 


-9 


Mazz05 


SN 1994D b < c 


CN 


Probable 


-11 


P96 


SN 1999gp 


SS 


No 


-4.5 


M08 


SN 1995D 


CN 


Uncertain 


+0 


S96 


SN 2000E 


SS 


Possible 


-6 


V03 


SN 1996X b < c 


CN 


Probable 


-1 


S01 


SN 2001V b 


SS 


Possible 


-7 


M08 


SN 1997br b 


SS 


No 


-9 


L99 


SN 2001el 


CN 


Probable 


-9 


Matt05 


SN 1997dt 


CN 


Uncertain 


-10 


M08 


SN 2003cg 


CN 


Possible 


-8.5 


E-R06 


SN 1998V 


CN 


Uncertain 


+0.5 


M08 


SN 2003du b ' c 


CN 


Probable 


-11 


S07 


SN 1998ab 


SS 


No 


-7.5 


M08 


SN 2005cf 


CN 


Uncertain 


-12 


W09 


SN 1998aq b 


CN 


Probable 


-9 


B03 


SN 2005cg 


SS 


No 


-10 


Q06 


SN 1998bu 


CN 


Probable 


-7 


H00 


SN 2005hj 


SS 


No 


-6 


Q07 


SN 1998es 


SS 


No 


-10 


M08 


SN 2006D 


CN 


Definite 


-7 


T07 


SN 1999aa 


SS 


No 


-11 


Gara04 


SN 2010ai 


CN 


Probable 


-10 


this work 


SN 1999ac b 


SS 


Possible 


-15 


Gara05 


PTFlOicb 


CN 


Probable 


-10 


this work 


HVG 


SN 1981B 


BL 


No 





B83 


SN 2000fa 


BL 


Uncertain 


-10 


M08 


SN 1984A 


BL 


No 


-7 


B89 


SN 2002bo 


BL 


No 


-13 


S05 


SN 1992A 


BL 


Probable 


-6.5 


K93 


SN 2002dj 


BL 


No 


-11 


P08 


SN 1997do 


BL 


No 


-11 


M08 


SN 2002er 


BL 


Possible 


-11 


K05 


SN 1998dh 


BL 


No 


-9 


M08 


SN 2004dt b 


BL 


No 


-7 


W06 


SN 1998ec 


BL 


No 


-2.5 


M08 


SN 2006X 


BL 


Possible 


-7 


Y09 


SN 1999cc 


BL 


No 


-3 


M08 


SN 2007gi 


BL 


Possible 


-7.5 


Z10 


SN 1999cl 


BL 


No 


-7.5 


M08 


SN 2009ig 


BL 


Possible 


-14 


Mil, in prep. 


SN 1999ej 


BL 


No 


-0.5 


M08 












FAINT 


SN 1986G 


CL 


Possible 


-6 


P87 


SN 2000cn 


CL 


Uncertain 


-9 


M08 


SN 1989B 


CL 


Possible 


-7 


W94 


SN 2000dk 


CL 


Possible 


-4.5 


M08 


SN 1991bg 


CL 


No 


+1 


T96 


SN 2004eo 


CL 


No 


-11 


P07 


SN 1997cn 


CL 


No 





T98 


SN 2005bl b 


CL 


Probable 


-6 


T08 


SN 1998bp 


CL 


No 


-2.5 


M08 


SN 2006bt 


CL 


Probable 


-3 


FlOb 


SN 1998de 


CL 


No 


-6.5 


M08 


SN 2010Y 


CL 


Definite 


-7 


this work 


SN 1999by b - c 


CL 


Possible 


-5 


Garn04 












SC d 


SN 2003fg 




Definite 




H06 


SN 2006gz b 




Definite 


-14 


H07 


SN 2007if 




Definite 


-9 


S10 


SN 2009dc 




Definite 


-8 


unpublished 


misc. 


SN 2000cx 




Possible 


-3 


L01 


SN 2005hk b ' c 




Possible 


-5 


S08 


SN 2002bj 




Definite 


+7 


P10 


SN 2007qd 




Possible 


+3 


M10 


SN 2002cx 




Probable 


-4 


L03 


SN 2008ha 




Definite 


-1 


F09, FlOa 



a SN type notation of iBranch et al . (2006), where CN = 'core normal', SS = 'shallow silicon', CL = 'cool', and BL = 'broad 
line'. The similar subtypes of lBenetti et al.l (|2005f ) have been used to separate to these SNe la into subclasses. 
b Taken from The Supernova Spectrum Archive, SuSpect. 

c Objects wi t h bot h publish ed spectropolarimetr y data an d C II A6580 ab s orptio n signatu res: 1994D;[W arig ct alj l|1996f ). 1996X; 
IWang et all (fl997h . 1999bv; IHowell et all IpOOH) . 2003du; ILeonard et all IpOOH ). 2005hk; IChornock et al.l |2006T ). 
d SC = Super-Chan dra Candidate s 

e Ref erences- (B89) IBarbon et al.l [19891: (B83) IBranch et aT]fT98o1; (B03) IBranch ct al. 2001 (E-R0 6) lElias-Rosa eFaLI [20061: 
(F99 ) iFisher et al.lll999l: (F09 ) Folcv ct al. 2009; (FlOa) Folcv ct al. 2010a|; ( FlOb) iFolevet al. 2010b; (Gara04) Garavin ilFaLl 
12001 (Gara05) IGaravini et al]|200l; (Garn04) IGarnayich et all 12001 (HOP ) Hcrnand eTeTaLl 120 00 : (H07) Hickcn etldl 120071 ; 
(H06) IHowell e t al. 2006;; ( K93) iKrrshner et al l fl99l : (K05) iKotak et all [20051 : (L99) iLi et all fl999l: (L01 ) iLi et all l200ll: 
ILi et all 120031; (Mil) IMarion eFaLI [20TH; (M08) iMatheson et alJ 120081 : (Matt05) IMattila et all [20051: (M01) IMazzall 
_, (Mazz05) iMazzali et all 120051: (M10) iMcClelland et al] [20101; (Mil) Mario n et al. 2011; (P07) iPastorello et alJ 120071: 
(P96) Patat ct al. 1996; (P87) Phillips ct al. 1987; (P08) Pignata ct al. 2008; (P10) Poznanski ct al. 2010; (Q06) Quimby ct al 



(LOS) 
[2001 : 



2006; (Q07) Quimby ct al. 2007; (S08) Sahu ct al. 2008; (S01) Salvo ct al. 2001; (S96) Sadakanc ct al. 1996; (S10) Scalzo ct al 



20101; (S07 ) Stanishcv ct al. 2007; (S05) Stchlc ct al. 2005 j; (T08) [T aubcnbcrger ct al. 200 8j; (T07) [ Thomas ct al. 2007; (T96 
Turatto et al. 1996; (T98) iTuratto et a l. 1998; (V03) Val entini et"aI1 120031; (W0 6) IWang et al.ll2006l : (W09) I Wang et al.lf200T 
(W94) Wells ct al. 1994; (Y09) [Yamanaka ct al. 20091; (ZlO) IZhang et al.|[20lol . 
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Table 3 

SYNOW Fit Parameters 



Supernova 


Epoch 


Vphot 


C II a 




Si II a 




Name 


(days) 


(km s —1 ) 


Velocity 


T 


Velocity 


T 


bJN 199U1N 


-14 


13,000 


16,000 


0.60 


>14,U0U 


0.50 




-8 


11,000 


14,000 


0.25 


11,000 


1.40 




-2 


11,000 


14,000 


0.15 


11,000 


16.0 




2 


11,000 


14,000 


0.35 


11,000 


16.0 




-11 


14,000 


14,000 


1.30 


>14,UUU 


14.0 




-10 


13,000 


14,000 


0.45 


>lo,UUU 


12.0 




-8 


13,000 


14,000 


0.45 


>13,000 


6.00 




-5 


12,000 


14,000 


0.20 


— i o r\r\r\ 

>12,000 


6.00 


bJN 199oX 


-4 


13,500 


13,500 


0.60 


13,500 


15.0 




-2 


13,000 


13,000 


0.40 


13,000 


13.0 




-1 


12,000 


13,000 


0.04 


13,000 


8.00 


blN lyyoaq 


-9 


12,000 


15,000 


0.20 


13,000 


2.10 




-8 


12,000 


14,000 


0.20 


13,000 


2.10 







12,000 


13,500 


0.60 


12,000 


9.00 


blN lyyyac 


-15 


13,000 


17,000 


0.40 


15,000 


2.00 




-9 


14,000 


17,500 


0.30 


14,000 


3.00 


biN lyyyby 


-5 


12,000 


14,500 


0.50 


12,000 


12.0 




-4 


12,000 


14,500 


0.30 


12,000 


13.0 




-3 


12,000 


14,500 


0.30 


12,000 


14.0 


SN 2001V 


-7 


8000 


12,000 


0.15 


8000 


1.60 




-6 


8000 


12,000 


0.15 


8000 


1.60 




-4 


8000 


12,500 


0.15 


8000 


1.80 


blN zUUoau 


-13 


14,000 


14,000 


0.60 


> 14,000 


>4.00 




-11 


12,000 


14,000 


0.60 


>lz,U0U 


>4.50 


blN zUUbtlk 


-5 


6000 


7500 


0.07 


7000 


0.30 







6000 


7500 


0.14 


^UUU 


0.70 


bJN zUUoD 


-7 


11,000 


14,500 


0.30 


13,000 


0.70 




-5 


11,000 


14,500 


0.20 


13,000 


0.70 


blN zUUobt 


-4 


6500 


6500 


0.30 


> 11,000 


>0.50 




-3 


6500 


6500 


0.25 


>11,000 


>0.70 


bJN iUUbgz 


-14 


12,000 


18,500 


0.50 


16,500 


1.10 


-12 


12,000 


18,500 


0.40 


16,500 


1.50 




-10 


12,000 


17,500 


0.30 


16,000 


1.60 


SIN 200711 


-9 


8500 


9500 


0.25 


8500 


0.30 


cm rinnoi c 

bJN zUUona 


-1 


1000 


1000 


1.00 


2500 


10.0 


blN zUUyac 


-8 


8000 


11,000 


0.50 


10,000 


0.80 




-7 


8000 


10,000 


0.30 


9500 


0.50 




-6 


8000 


10,000 


0.30 


9500 


0.50 




-5 


7000 


9500 


0.30 


9000 


0.50 




-3 


7000 


9500 


0.30 


8500 


0.50 




-2 


6500 


9000 


0.30 


8000 


0.50 




-1 


6500 


9000 


0.30 


8000 


0.50 







6500 


8500 


0.40 


8500 


1.10 




2 


6500 


8500 


0.30 


8000 


0.60 




3 


6000 


7500 


0.30 


7500 


0.70 




4 


6000 


7500 


0.30 


(oOU 


0.70 




5 


6000 


7500 


0.30 


7500 


0.70 






5500 


7500 


0.30 


f UUU 


0.70 




7 


5500 


7500 


0.30 


6500 


0.70 




8 


6000 


7500 


0.30 


7000 


0.70 




Q 


ouuu 


uuuu 




7finn 


1 no 




15 


5000 


5000 


0.20 


6500 


0.80 


SN 2009ig 


-14 


16,000 


16,000 


1.00 


>16,000 


>24 




-13 


16,000 


16,000 


0.50 


>16,000 


>24 


SN 2010Y 


-7 


12,000 


15,000 


0.50 


14,000 


6.00 




-(5 


12,000 


15,000 


0.50 


14,000 


5.00 




-3 


12,000 


14,500 


0.50 


13,000 


10.0 




-2 


11,500 


14,000 


0.50 


12,500 


18.0 




1 


11,000 


14,000 


0.20 


12,000 


20.0 


SN 2010ai 


-10 


13,500 


14,000 


1.00 


13,500 


2.00 




-8 


12,500 


14,000 


0.50 


12,500 


1.50 


PTFlOicb 


-10 


11,000 


14,500 


0.10 


>12,000 


>1.10 



a Fits made with Tea^lO.OOO K and v e = 1000 km s . 

b Us ed two-components o f Si II to fit broad Si II A6355 P-Cygni profiles. 

c See lValenti et al. (2009) for an alternate SYNOW fit where the origin of explosion is interpreted to be a core-collapse supernova. 



